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Abstract 

The electric-double-layer structure in an electrolyte close to a solid substrate near the three-phase contact line is approximated 
by considering the linearized Poisson-Boltzmann equation in a wedge geometry. The mathematical approach complements 
the semi-analytical solutions reported in the literature by providing easily available characteristic information on the double 
layer structure. In particular, the model contains a length scale that quantifies the distance from the fluid-fluid interface over 
which this boundary influences the electric double layer. The analysis is based on an approximation for the equipotential lines. 
Excellent agreement between the model predictions and numerical results is achieved for a significant range of contact angles. 
The length scale quantifying the influence of the fluid-fluid interface is proportional to the Debye length and depends on the 
wall contact angle. It is shown that for contact angles approaching 90° there is a finite range of boundary influence. 
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The wetting of solid substrates by electrolyte solutions plays an important role in microti ui dies, specifically in digital 
microfluidics based on the electrostatic actuation of sessile droplets [T] [2]. In this context, it is important to study the 
influence of electric fields around the three-phase contact line. In electrowetting-on-dielectric, it has been shown that 
the large magnitude of the electric field close to the contact line can be the cause of contact-angle saturation |3J H] . 
Also without applying an external electric field various challenges remain, which have been addressed in Refs. [5HH] 
focusing on the electrostatic contribution to wetting. It has been shown that the energetic contribution of the electric 
double layer close to the three-phase contact region represents a significant fraction of the total free energy of a sessile 
droplet that is moderately sized compared to the Debye length |SHH]. However, little is known about the structure of 
the electric double layer in the contact-line region, a topic we address in the present study. Hence we consider the 
three-phase contact region of two immiscible fluids, e.g. water (index w) and oil (index np for nonpolar), and a solid 
substrate (index s); see Fig. la The water-substrate interface is assumed to carry a constant surface charge density p s 




Figure 1. Schematics of the physical problem and the approximation: (a) geometry, material properties, and boundary 
conditions for the potential <f> and | (b) | principle of approximation, parametrization of the functional form of the equipotential 
lines, and definition of the length scale C. 

corresponding to a dielectric substrate, whereas charges are supposed to be absent at both the oil-substrate and the 
oil-water interface. All interfaces are assumed to be infinitely thin and uncurved, constituting a wedge geometry for 
the water phase. We therefore investigate the system on a length scale large enough to neglect the effects of long-range 
van der Waals forces |10H12| . but small enough to neglect the interfacial curvature due to the electric field [5] or 
gravitational forces. In equilibrium, the electrostatic potential distribution cf> within the water phase (assumed to 

be a 1-1 electrolyte), can be modeled by the Poisson-Boltzmann equation A</> = kT/eXjj sinh (e<p/kT), where k is the 
Boltzmann constant, T is the absolute temperature, and e is the elementary charge. In the case that the electrostatic 
energy of the ions is much smaller than their thermal energy, implying that e<p/kT « 1, the Poisson-Boltzmann 
equation can be linearized to yield 



A0=t^ (1) 
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which is commonly referred to as the Debye-Hiickel approximation. For general values of the dielectric permittivities e w , 
e np and e s , the potential in all of the three phases would have to be calculated because of the coupling through the 
interfacial conditions. The corresponding conditions for the normal components of the electric field can be written 
as n • (e w V4> w - £i^4>i) = Pi, where i stands for s or np, n is the normal vector pointing out of the water phase and pi 
the surface charge density at the respective boundaries |13) . with p np = 0. If we assume that £i/s w « 1, i e [s,np] 
corresponding to the comparatively high permanent dipole moment of water molecules, the interfacial conditions 
are simplified to n-S7(j) w « Pi/e w (again i e [s,np]; cf. Fig. |la[ ). This simplification allows for a decoupling of the 
electrostatic potential in the water phase from the potentials in the remaining phases, so that we may focus solely 
on the water phase. At this point, we introduce nondimensional quantities by scal ing the potential (f> with p s Xd/sv 



and the coordinates x, y and f with the Debye length Ad, respectively ([cf. Fig. la ]. Then Eq. (JlJ and the boundary 
conditions for the problem transform into 

A0 = «£ (2) 

0Uoo = (3) 

Vf4U oo = (4) 

V(f>-e e \ e=0 = -l (5) 

V0-e e | e=Q+7r/2 = O (6) 
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Eq. Q together with the boundary conditions ([3])-([6]) can be solved analytically using the Kontorovich-Lebedev 
transform [5HH1 114) 115] . This method provides representations of the solution in the form of one-dimensional integrals 
that are usually evaluable only by means of costly numerical quadrature. In general, such integral representations do 
not reveal the functional dependence of the solution on physical parameters, for example, allowing for the study of the 
decay behavior or the uncovering of characteristic length scales. Only an asymptotic evaluation of the integrals is 
possible |5j [14] , yielding approximate solutions of limited validity. 

In this study, we develop a semi heuristic approximation of the electrostatic potential distribution described by 
Eq. (Q in a three-phase contact region. The paper is organized as follows. First an approximation for the shape of 
the equipotential lines is developed and compared with numerical results. From this model, a characteristic length 
scale over which the oil-water interface influences the double layer structure is extracted. Since the derived expression 
is implicit with respect to the potential, the latter needs to be evaluated by means of Newton's algorithm using a 
procedure described at the end of the paper. 

Before developing the approximation, it is instructive to analyze the general properties of the solution as they can 
be readily extracted from the boundary conditions ([3])-([6]). Condition (|3| forces the potential to vanish at y -> oo, 
whereas ([5| determines the normal electric field by prescribing the surface charge density. Condition (JsJ) implies 
that the electric field vector is parallel to the oil-water interface for 6 = a + 7r/2, whereas condition Q expresses 
the fact the electric field becomes normal to the substrate surface for x -*■ oo, approaching asymptotically the far 
field solution cj)(x -> oo,y) = exp(-y). In addition, the Green's function of Eq. (JlJ has a dominant exponential decay 
characteristic [T3]- Therefore, regarding the influence of the contact line region as a perturbation to the electric double 
layer extending into the half-space with positive values of x, it is reasonable to assume an exponential transition between 
the potential in the far-field and close to the contact line that occurs over a characteristic length scale describing 
the range of the boundary influence. This behavior can be captured directly by approximating the shape of the 
equipotential lines through a function y = f(x). In order to parametrize f(x), we introduce a pair of coordinates (770, yo) 



as depicted in Fig. lb where yo is measured at x -*■ 00. Note that f(x) also depends on a as well as on yo though we 



do not explicitly include this dependence in the notation for brevity. For the function /(x) we choose the ansatz 

f(x) = A-Be- Cx (7) 

subject to the conditions 

df 



dx 



= tana (8) 

x=—rio sin a 

f{x - 00) = y (9) 

f{x = -rjo sin a) = 770 cos a (10) 

where condition (|8| corresponds to ^ and condition (|9| corresponds to Q, respectively. Combining the ansatz ^ 
with the conditions ([8]) ( 10 1 yields 



f(x) = ?/o -(2/o-?7o cos a) e '"°»»-™-° (11) 
Since the choice of the potential value completely determines y for each equipotential line through the far-field relation 

(j) = e- yo (12) 



it only remains to find a relation between y and 770 in order to eliminate 770 from Eq. (111. For this purpose we utilize 
a numerical solution of Eq. Q subject to the boundary conditions (|3j — (J6j) obtained by means of the commercial finite 
element solver COMSOL Multiphysics®. The parameters of the numerical calculations, especially the grid structure 
and size, were chosen such that grid-independent results were obtained. Fig. [2] shows the results for different values 
of a. Obviously, a linear relation yo Vo m &y serve as an excellent approximation for -0.2-7T < a < 0.2-7T, whereas 
linearity is slightly violated for |a| > 0.27T. Therefore we assume that 

y = ci] + d (13) 



where the dependence of c and d on a is suppressed in the notation. Eq. (13) in conjunction with Eq. (12) implies 



that the potential distribution at the oil-water interface tfi ow can be approximated by the function 

4>*w = e- cr >°~ d =: 6e- c * (14) 

In the next step, we propose expressions for the coefficients c and d (or b) from analytical considerations and compare 
them to numerical results. The relation for b is found from an asymptotic evaluation of the integral representation of 
the potential (j) at the point (x,y) = (0,0) to be of the form |5j|6] 

cfiowivo = 0) = <K* = 0, V = 0) = b = e- d = (15) 

7r + 2a 
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Figure 2. Relation between the coordinates yo and 770 for different values of a extracted from numerical results 



Furthermore, the gradient of <fi a t ( x ,y) = (0,0) is parallel to the oil-water interface [cf. Fig. lb ] because of the 
boundary condition ([6]), implying that its absolute value is given by \d(f> ow / dr]o\. In addition, the gradient has to fulfill 
condition |5| . Therefore it follows from a comparison between the boundary conditions for (j) at x -*■ and 00 that 



1 



cos a ch/o 



yo=o 



dr] 



1 



£=0, ?7o=0 



7r + 2a 



b cos a 7r cos a 



(16) 



In order to compare the expressions (H5| and (16 I with the numerical solution, we fit Eq. (14 1 to the numerically 



calculated potential distribution at the oil-water interface and plot the resulting coefficients b and c, as well as their 
proposed dependencies on a, in Fig. [3} While Eq. (15) predicts the values accurately [Fig. 3a ], expression (16 1 turns 
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Figure 3. Comparison of the forms ( 15 1 and ( 16 \ for b and c with coefficients found from fitting Eq. ( 14 1 to the numerical 
solution 



out to be inappropriate for modeling the coefficient c [Fig. 3b ], which can be reproduced by the linear fit c = 1 + 0.36a 
instead. A careful investigation reveals that, although Eq. (16) correctly predicts the gradient of (j) at (x,y) = (0,0), 



the intention to construct an accurate approximation to the entire potential distribution at the oil-water interface 
rather demands a linear relation between c and a. This contradiction is caused by the fact that, except for a = 0, the 
potential distribution at the oil- water interface deviates from an exponential function. As a result, the function / 
modeling the shape of the equipotential lines is given by 



r/ \ / V0~d \ -tana 

j{ x ) = yo-\yo cosale =«o-(yo- d ) = 

(7r + 2a\ 
j ; c = 1 + 0.36 a ; yo = - ln< 



(17) 



Since up to this point we have only made sure that the model reproduces the potential at the oil-water interface, it is 
necessary to compare the equipotential lines in the full domain to the numerical results. Examples are depicted in Fig. [4] 
for four different values of a. The agreement within the range \a\ < 0.25n/2 is excellent, confirming that the shape of the 
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(a) a = -0.25tt/2 



(b) a = -0.17I-/2 
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Figure 4. Comparison of the equipotential lines (solid lines) calculated from Eq. 1 17 1 to the numerical results (symbols) for 
several values of a; lines and data points, from bottom to top, correspond to <j> = 7r(10/ll, 9/11, . . . , l/ll)/(7r + 2a) for a < and 
to <f> = (10/11, 9/11, . . . , 1/11) for a > 0, respectively 



equipotential lines in this range can be approximated by exponential functions. For larger angles 0.257r/2 < \a\ < OAtt/2 
(not shown in this paper), Eq. (17) still provides a good approximation to the equipotential lines exhibiting only small 
deviations. Beyond this region, the agreement between model and numerical results deteriorates with increasing values 
of \a\, but at least qualitative agreement can be observed for OAir/2 < \a\ < 0.57r/2. Note that the linear behavior y oc ?y 
shown in Fig. [2] holds for values up to |a| < 0.47r/2 at least 



As a main result, from Eqs. (Ill and (17 1 we can extract a characteristic length scale that is given by [cf. Fig. lb 



yo~VoCOsa j/o ( c ~ cos a) + dcosa 
C = Ad = - \d 



tana 



ctana 



(18) 



The length scale C in Eq. (18) is proportional to the Debye length Ad, but is also dependent on the angle a and the 



far- field coordinate yo , and describes how far the influence of the contact- line region extends into the water phase until 
the electric double layer finally reaches its far-field equilibrium structure. For small values of a we find 



A 



D 



+ 0.36 y + [" A (1 + 0.36 tt) + ( - - 0.36 2 ) 
. it V 2 / 



a + 0(a 2 ) 



(19) 



Expansion ( 19 1 reveals an important property of the potential field, namely a finite value of the length scale C appearing 
even for a -* 0. While for vanishing a the prefactor of the exponential term in Eq. (17 1 also vanishes and thus cancels 



the effect of £, there is always an exponential decay over a finite length scale for arbitrarily small a ^ 0. 

Besides the characteristic length scale C in Eq. (18 1, the approximation (17) also yields an implicit expression of 
the form y = /(a;,0), which cannot be solved for </> analytically. Therefore, if the potential at a certain point shall be 
evaluated, numerical methods are needed. Note that in contrast to the complex numerical evaluation of an integral 
representation of the analytical solution [S], the numerical method to solve Eq. (17) may be the inexpensive Newton 
algorithm for nonlinear equations |16| . The evaluation of the potential </> at a given point (x,y) requires the following 
steps. Therein, we exploit the much smaller variation of the field y (x,y) with the coordinate y compared to (j)(x,y), 
whose exponential variation is removed by Eq. (12). 
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1. Estimate a starting value for yo using the far- field relation yo = V as x 00 



2. Solve Eq. ( 17) for yo iteratively by means of Newton's algorithm which yields fast convergence within only a few 
iterations 



3. Calculate <f>{x,y) from Eq. (12) 



In this way, the model ( 17 1 not only provides characteristic information on the double layer structure but also allows for 
an inexpensive calculation of the complete potential distribution. Note that by using the nonlinear Poisson-Boltzmann 
equation instead of the Debye-Hiickel approximation a similar model can be obtained, though the value of the 
electrostatic potential on the three-phase contact line is not known in this case. 
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